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Using Kepler transit observations to measure stellar spot belt 
migration rates 
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ABSTRACT 

Planetary transits provide a unique opportunity to investigate the surface distributions of star 
spots. Our aim is to determine if, with continuous observation (such as the data that will be 
provided by the Kepler mission), we can in addition measure the rate of drift of the spot 
belts. We begin by simulating magnetic cycles suitable for the Sun and more active stars, 
incorporating both flux emergence and surface transport. This provides the radial magnetic 
field distribution on the stellar surface as a function of time. We then model the transit of a 
planet whose orbital axis is misaligned with the stellar rotation axis. Such a planet could occult 
spots at a range of latitudes. This allows us to complete the forward modelling of the shape of 
the transit lightcurve. We then attempt the inverse problem of recovering spot locations from 
the transit alone. From this we determine if transit lightcurves can be used to measure spot 
belt locations as a function of time. We find that for low-activity stars such as the Sun, the 3.5 
year Kepler window is insufficient to determine this drift rate. For more active stars, it may be 
difficult to distinguish subtle diff'erences in the nature of flux emergence, such as the degree of 
overlap of the "butterfly wings". The rate and direction of drift of the spot belts can however 
be determined for these stars. This would provide a critical test of dynamo theory. 

Key words: Stars: Activity, Stars: Spots, Stars: Planetary Systems; Planet and Satellites: 
General; Sun: Sunspots. 



1 INTRODUCTION 

The solar magnetic cycle is now very well observed. On the Sun, 
the distribution of solar spots is known to vary over an 11 -year cy- 
cle. The latitudes at which spots emerge as a function of time gives 
rise to the well-known "butterfly diagram". This diagram shows 
that at the beginning of a solar cycle the emergence of spots prefer- 
entially occurs at high latitudes and moves down to lower latitudes 
towards the end of the cycle. It was initially assumed that this would 
be the case on all stars, however, stars such as AB Doradus, a very 
bright, rapid-rotator have be en observed over many years wit h no 
such pattern being detected jPonati & Collier Cameron|[l997h . By 
monitoring the variability in stellar brightness over time we can 
build up a picture of the spot dis tributions th a t caus e rotational 
modulations in stellar lightcurves. iBerdvuginal (l2005h provides a 
comprehensive review on star spot theory and observations on cool 
stars. 

The transit method of detecting exoplanets has now advanced 
well beyond the initial detection of hot-Jupiters, with observations 
revealing more information about both the planet and the host star. 
Asymmetries in exoplanetary transit lightcurves are proving to be 
a useful tool in probing both the planetary atm osphere and the in - 
teraction between the planet and its host star. IFossati et 
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reported an asymmetr y in the light curve of WASP- 12b when ob - 
served in the near-UV. Ivrdotto et al.l ilOld) and lLlama etaP (1201 ih 
have shown that this asymmetry could be caused by a magnetically- 
channelled stellar wind interacting with the planet's magnetic field. 
This results in a bow shock transiting ahead of the planet, suggest- 
ing that transit observations may be an indirect method for detect- 
ing exoplanetary magnetic fields. 

Another reported asymmetry in transit lightcurves arises when 
a planet transits over a star spot, as shown in Figure [T] When the 
planet is transiting over a dark region on the stellar disk such as a 
star spot, the amount of light blocked by the planet be comes lower 
which produces a positive "bump" in the lightcurve jRabus et al.l 
[2009,) . There have been reports of planets transiting across star 



spots i n many different systems (see for example [Rabus et al 



2009| ) : IPont etaPfeOOThilLanza et al.l ( [200^ : ISanchis-Oieda et al. 
201 lh ). ISilva- Valid feOOsI) have shown that detecting star spots 



through transit observations allows us to determine the longitude 
of the spot and therefor e multiple transit events can help determine 
the stellar rotation rate. ISanchis-Oieda et aP (l201lh have also used 
such transits to investigate the alignment between the stellar rota- 
tion axis and the planet's orbital axis by looking for recurring star 
spot bumps in the transit lightcurves. Whilst most of the work so 
far has concentrated on individual spots, missions such as Kepler 
provide a unique opportunity to build up a picture of spot evolu- 
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tion and cycles on planet hosting stars by providing highly precise, 
continuous observations of many transiting systems. 

In this Letter we show that it may be possible to measure the 
migration rates of stellar activity belts from transit observations 
such as those provided by Kepler. Using simulated solar and stellar 
cycles we investigate whether this method can distinguish between 
different types of magnetic cycle. 

2 THE MODEL 

2.1 Simulated butterfly diagrams 

In order to study how transit light curves may help recover the spot 
evolution on planet-hosting stars, we first specify the emergence of 
flux through the stellar surface. The choice of the emergence rates 
and strengths of the bipoles, as well as the time-latitude dependence 
are i mportant elements to accurately simulate the evolution of star 
spots (iMackav et al.ll2004l) . 

For this investigation we choose various examples of possible 
solar and stellar butterfly diagrams. Table [T] provides an overview 
of our four models and the left column of Figure[2]shows the corre- 
sponding input fl ux emergence patterns. Mode l A is the solar but- 
terfly diagram of Ivan Ballegooiien et al.l (Il998h . As shown in Fig- 
ure[2al this pattern of magnetic flux emergence reproduces the solar 
11 -year spot cycle. The beginning of the cycle shows the bipoles 
emerging preferentially at latitudes of +40°. The bipoles emerge at 
lower latitudes through the cycle, finally reaching latitudes of + 10° 
by the end of the cycle. 

Models B, C, and D ar e based on theoretical butterfly dia- 
grams used by iMcIvor et al.l (l2006h to simulate magnetic cycles 
on active stars such as AB Dor. AB Dor is a solar-like rapid lo- 
tator (Prot - 0-51 days), with a stron ger magnetic field than the 
Sun toonati & Collier Cameroniri997h . Unlike the Sun, which ex- 
hibits clear belts of activity whose latitude varies through the mag- 
netic c ycle, AB Dor show s spots at all latitudes at each observing 
epoch irefl'ers et al]|2007h . The flux emergence patterns used by 
I all ( l2006h ir 



I in conjunction with a model for the transport 



iMcIvor et 

of this flux across the stellar surface reprodu ce this o bserved be- 
haviour. For ease of comparison with the Sun jMcIvoT et al. ( 2006,) 
assume an 11 -year cycle in their simulated butterfly diagrams, 
which will also be adopted here. Model B, shown in Figure [2cl 
consists of an enhanced butterfly pattern, with the maximum lat- 
itude of emergence extended to 70°. Model C (Figure [Je]), is an 
overlapping butterfly pattern, again with the maximum latitude of 
emergence increased. This input allows flux to appear at both high 
and low latitudes at all times but also retains the butterfly pattern. 
Model D (Figure[2g]), assumes no latitude-time relationship and the 
emerging bipoles are completely random in both longitude and lat- 
itude throughout the cycle. 

2.2 Surface evolution 

Once the magnetic flux has emerged through the stellar surface it 
will be transported across the stellar surface by difl'erential rota- 
tion, meridional flow and difl'usion. It is the coupling of the flux 
emergence with these transport processes that results in a cyclic 
behaviour of the surface flux. To reproduce the solar and stellar 
cycles we use a mag netic flux transport code as described fully in 
iMackav et al.l (|2004|) (and references therein). By accounting for 
the efl'ects of difl'erential rotation, meridional flow and super gran- 
ular difl'usion, we evolve the stellar magnetic field at the surface of 
the star, i.e. r - R^. The model only accounts for the large-scale 
magnetic field with spatial scales > 30Mm, but has been shown 
to accurately reproduce the solar activity cycle (see for example. 




Figure 1. Schematic of the model (not drawn to scale). The planet transits 
the star on an inclined orbit. The angle (3 defines the obliquity of the planet. 
When the planet transits over a dark region of the stellar disk such as a star 
spot, a bump appears in the transit lightcurve. 



|Wang & Sheeley (1991)). The simulated butterfly diagrams (Fig- 
ures[2al|2c][2e]and[2g]) are taken as input to the code which evolves 
the radial magnetic field Br(Ri,, 6, 0, by solving the partial diff'er- 
ential equation: 



dBr 



1 d 

sinOde 



sin 6 



dBr 

-u{e)Br + D-^ 
do 



D 



where D is the photospheric diff'usion constant, Q.{6) is the profile 
for differential rotation and u{6) describes the meridional flow. The 
parameters used for each model are given in Table [T] 

No direct measurements of meridional flow have been made 
on main sequence stars. We therefore assume a solar profile, where 
the flow is directed poleward and is expressed in terms of latitude 
{X = nil - 6)hy the function: 



u{X) ■■ 



-uq sm(7TA/Ao), 



0, 



otherwise. 



(1) 



Under this prescription, the meridional flow velocity vanishes 
above Aq. We use th e observed values for the Sun of Ap = 75° 
and Wo = 11 ms"^ (ISnodgrass & Dailevlll996l : iHathawavl Il996h 
for the solar type s tar and uq = 100 ms~^ for the rapid rotator 
(iMackav et al.ll2004l) . For all models, the difl'erential rotation is ex- 
pressed using the solar profile of (Snodgrass 1983), 



Q(e) = 13.38 - 2.30 cos^ 6 - 1.62 cos^ 6 - Qq deg day" 



(2) 



where Qq is the Carrington rate (13.20 deg d"^. Finally, the super 
granular diff'usion constant D, determines the timescale on which 
flux moves from the centre of granular cells to the boundaries where 
it can interact and cancel with flux of opp osite polarity. T his is 
given the solar value of Z) = 450 km^ s~^ (lLeightonlll964l) . The 
evolution code outputs a radial magnetic map at each time step, 
which is set here to be three days to simulate the typical times be- 
tween a transit event. 

2.3 Simulated lightcurves 

To produce the lightcurves we assume a star of radius = Rq and 
a hot- Jupiter of radius Rp = Rj, so that a 1% dip in flux will be 
recorded as the planet transits across the stellar disk (assuming the 
portion of the stellar disk occulted by the planet has intensity 0). 
A schematic of the setup is shown in Figure [T] We take the impact 
parameter, b - so that the planet transits across the centre of 
the stellar disk. The angle between the stellar rotation axis and the 
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Figure 2. The left column shows the input flux emergence patterns for each of our models (Ivan Balleeooiien et al.ll998l : lMcIvor et al .l20Q6h . The right column 
shows the distribution of spots recovered through the transit method as explained in Section 2.3. The top row corresponds to Model A, a Solar type star. The 
final three rows correspond to Models B, C, and D which represent the rapid-rotator with an enhanced butterfly pattern, an overlapping butterfly pattern and, 
no butterfly pattern respectively. 
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Table 1. Parameters used for each of our models. For each model Ri, = Rq, 
Rp = Rj, and /? = 30° (see Figure [TJ. The second column gives the type 
of input butterfly pattern used; Column three is the initial meridional flow 
speed; Column four is the diff'usion constant; Column five is the drift rate 
of the input latitudes of flux emergence; Column six is the drift rate of the 
predicted spot latitudes. 

Model Type uq D P^ot Input drift Recovered drift 

(ms~^) (km^ s~^) (days) (degyr~^) (degyr~^) 

A Solar 11 450 27 -2.2 -2.0 
butterfly 

B Rapid rotator 100 450 0.5 -4.8 -4.5 
butterfly 

C Rapid rotator 100 450 0.5 -3.3 -2.5 
overlapped 

D Rapid rotator 100 450 0.5 
no butterfly 



orbital plane of the planet (as shown in Figure [B is taken to be, 
- 30°. Under this setup, the planet transits over a large range of 
latitudes on the stellar disk and also spends an equal amount of time 
in each hemisphere. If the planet were to spend more time transiting 
over one hemisphere of the stellar disk an observational bias may 
be present in the recovered distribution of spots. We also assume 
that Porb = 3 days. We set the stellar rotation period to be Prot = 27 
days for the solar type star and Prot = 0.5 days for the rapid rotator. 
In all cases the star is edge on to the observer. 

At every transit crossing event (i.e. once every three days) our 
surface flux transport simulations output the magnetic configura- 
tion on the surface of the star allowing for both the rotation of 
the star and also foreshortening. From this we classify regions of 
strongest magnetic field (|B| > 0.9 IBIj^ax) as spots. 

We irradiate the stellar disk with Monte Carlo photon packets 
to reproduce the spatial intensity distribution of a limb-darkene d 
star. We use the limb-darkening law as given bv IClared Ew^ , 
where the intensity / is given by 

M = l-f^^a.il-,"'\ (3) 

where ji = cos^ = (1 - r^Y^^, < r < 1. The coeflficients are 
chose n to match the limb darkening of Kepler observations (ISingI 
For this analysis, we choose the limb-darkening coefl&cients 
for a star with Teff = 7500 K and logg = 4.5. The spatial intensity 
distribution is also aflfected by the spot distribution. For simplic- 
ity, we assume that regions that contain a spot have an intensity 
of 0, i.e. are completely dark. We then sum the total radiated flux 
and normalise by the total number of photons to create the transit 
lightcurve. 

Any instances where the planet occults a spot, or part of a spot 
then the fractional loss of light will be less and a positive bump will 
be registered in the lightcurve. Should the planet occult multiple 
spots in the same transit we are able to recover each spot crossing 
separately. The phase of any bumps can then be converted back 
into latitude on the star. The presence of noise in the data of course 
li niits the sizes of spots that can b e recovered. Guided by the results 
of ISanchis-Oieda & WinnI (I20T 1|) we only recover spots where the 
deviation in normalised flux is larger than 0.001%. The phase of 
the bump provides the position of the planet (x,y) on the projected 
stellar disk. This location is then converted into latitude on the star 
using using the equation: 



^ lysini Jx^ +yA 
A = sm-^ ^ ^ . (4) 

By plotting spot latitudes as a function of time, we can therefore 
determine the rate and direction of drift of the spot belts. 

3 RESULTS 

The results of our simulations are shown in the right hand column 
of Figure [3 The recovered spot locations for Model A, a Solar 
type cycle are shown in Figure [2bl The results indeed show a time- 
latitude relation similar to the input flux emergence diagram. The 
recovered spots are at higher latitudes at the beginning of the cycle 
than at its end, with the pattern repeating after 1 1 -years. 

The recovered spots of the rapid rotator are just as revealing. 
Figure l2d1 shows the results for Model B, the enhanced butterfly 
diagram (without an overlap). The recovered spots show a cycle 
repeating on an 11 -year scale which matches the input. The spots 
appear at high latitudes at the beginning of the cycle and lower 
latitude spots appear towards the end of the cycle. Although the 
input pattern is reproduced, we notice high-latitude spots recov- 
ered throughout the cycle. This is to be expected due to the in- 
creased rotation rate of the star and also the increased meridional 
flow (co mpared to the solar values) that drives the magnetic field 
poleward (iMcIvor et al Even in the phase of the cycle where 

spots emerge at low latitude, they are rapidly pushed to higher 
latitudes because of the meriodional flow, which results in high- 
latitude spots being more readily recovered throughout the cycle. 

Figure [23 show the resultant spot latitudes for the overlapped 
butterfly diagram (Model C). This again shows that it is possible 
to recover a time-latitude dependence but suggests that it may be 
difficult to distinguish it from Model B (Figure l2dll. 

Finally, Figure |2h] shows the recovered spot latitudes for a 
cycle with no time-latitude dependence (Model D). Although the 
input diagram has no time-latitude dependence, the results show 
a weak latitude dependence, with more spots appearing at higher 
latitudes. This is likely due to the increased Coriolis force and 
meridional flow pushing the spots poleward. One interesting result 
is that the recovere d spots do not appear above 60° even though 
jMcIvor et al.ll2006b report magnetic activity at latitudes up to the 
pole. 

For this method to recover spots at very high latitudes, the 
planet must be on an highly inclined orbit (i.e. the value of /3 in 
Figure [T] must be very small). Foreshortening of the spots on the 
stellar disk causes high latitude spots to appear smaller than sim- 
ilar spots at lower latitude. This method is therefore less sensitive 
to detecting spots at high latitude and preferentially finds spots at 
lower latitudes than spots near the pole. 

To further compare the results of our simulations with the in- 
put, we calculate the drift rate of both the input spot emergence 
pattern through the cycle and the drift rate of the recovered spots. 
For the input flux emergence we are able to calculate the drift rate 
directly. However, since the recovered spots have been subjected to 
the surface evolution processes described in Section 2.2, the drift 
rate is recovered by measuring the slope of the lower envelope of 
each butterfly wing. This is achieved by taking 0.1 -year time bins, 
and then calculating the lower W'^ percentile of each bin. We then 
fit the function A(t) = A + Bt, where A is the initial latitude and 
B is the drift rate. The recovered value of the drift rates are shown 
in the final two columns of Table [T] In all four models the values 
compare favourably, however, the recovered drift rate is lower than 
the input rate. This is likely due to the meridional flow moving flux 
poleward. 
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Figure 3. Histograms comparing the drift rates of spot belts calculated from (a) the full simulation (red-line) with (b) the values that can be recovered from 
taking 10,000 3.5-year "snapshots" of the data which simulates the current Kepler observing window. Table [T] summarises the parameters used in each model. 



4 DISCUSSION 

In this Letter we have addressed the possibility of recovering spot 
cycles on stars hosting a transiting planet. This study is particularly 
timely because of the availability of the Kepler data. By modelling 
the emergence and transport of flux on the surface of a star through- 
out its magnetic cycle, we have predicted the distribution of spots 
as a function of time. From this forward-modelling we have demon- 
strated that transit lightcurves can be used to recover the drift rates 
of spot latitudes on active stars. This rate of drift of spot latitudes is 
a critical prediction of dynamo theories, but it has in the past been 
particularly difficult to observe directly. While other methods such 
as Doppler imaging can provide a snapshot of the spot distribu- 
tion, only continuous viewing capability such as that of Kepler can 
provide a sufficiently long time series to measure drift rates. The 
recovered drift rates of the active latitudes reveal the nature of the 
dynamo - whether it is a solar-like behaviour where the active lat- 
itudes drift towards the equator over the cycle, or anti- solar where 
they drift towards the poles. 

We have used various simulated spot cycles, assuming a cycle 
length of 11 -years, similar to that of the Sun. Kepler is currently 
scheduled to collect data for 3.5 years. To investigate whether it is 
possible to detect a drift in the recovered spot positions over this 
time we have randomly placed 3.5 year windows over each of our 
result sets. We then carry out the analysis as described in the final 
paragraph of Section 3 to determine the drift rate for each window. 
By repeating this 10,000 times we have investigated whether a con- 
sistent value for the drift rate can be recovered. 

Figure [3] shows the results for the data in the Northern hemi- 
sphere. We have carried out this analysis on Models A, B, and C 
(top three rows of Figure [2l). Model D (the bottom row of Figure 
|2]) is excluded from this analysis because it was designed to not 
exhibit any time-latitude dependence. 

For the Solar cycle (Model A) we find that a 3. 5 -year window 
of observation is not long enough to reveal a time-latitude depen- 
dence in the data. The solar case is hampered by the relatively low 
activity rate of the star, and hence, the likelihood of a spot being 
discovered through transit methods is significantly lower than for a 
more active star. To determine a time-latitude dependence for a so- 
lar type star a longer set of observations would likely be required. 

For more active stellar cycles (Models B, and C) the 3.5-year 
window produces a stronger Gaussian shaped distribution suggest- 
ing it is possible to recover a time-latitude dependence in the data. 
We find however that even if the general trend of spot locations and 
evolution can be recovered, it may be difficult to distinguish be- 
tween an enhanced spot cycle and an enhanced but overlapped spot 
cycle. 



The method of using transiting planets to determine changes 
in spot latitudes is particularly powerful, but requires not only that a 
planet is transiting, but also on an inclined orbit. Low latitude spots 
are preferentially recovered by this method, and shorter magnetic 
cycles will be easier to study within the Kepler lifetime. 

We have shown that the data collected using the observational 
window of the Kepler space mission has the potential to determine 
the drift rates of spot belts on active stars. Transit observations can 
therefore provide new insight into stellar activity and cycles, and 
hence are an additional, critical test for stellar dynamo theories. 
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